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Abstract. Over 800 sq. deg. of high Galactic latitude sky have been mapped at 21 
cm with the Robert C. Byrd Green Bank Telescope (GBT). An improved knowledge 
of the telescope's beam characteristics has allowed us to reliably map not only regions 
of high column density, but also such regions as ELAIS Nl, a targeted Spitzer field, 
which have very low HI column density. The additional fields we have observed cover 
a cross-section of dynamically and chemically interesting regions as indicated by the 
presence of intermediate/high velocity gas and/or anomalous far-IR (dust) colour. 



1. Introduction 



At high Galactic latitudes, the column density of ato mic hydrogen is known to correlate 
very strongly with far infrared emission from dust (IBoulanger et al.lll996l) . However, 

it has also been observed that HI clouds that have large peculiar velociti es (IVCs and 

HVCs) can have a significantly different FIR/HI ratio than normal disk gas ( Wakk er & Boulanger 
19861 : iHerbstmeier et al.lll993l : iMiville-Deschenes & Lagachdl2005l : iPeek et al.ll2009h . 
Deviations from the standard FIR/HI ratio are also sometimes found in isolated ISM 
components. Progress in understanding of the dust at high Galactic latitudes, and its 
connection to processes in the interstellar medium, hinges on an understanding of the 
gas-dust correlations at the various observed infrared and sub-mm wavelengths. With 
the Planck and Herschel satellites now returning sub-mm wavelength data at consid- 
erably higher angular resolution than the currently available all-sky HI surveys (the 
Leiden/ Argentine/Bonn (LAB) HI survey (Kalbe rla et al.ll2005 )). it is vital to have com- 
parable resolution, high sensitivity 21cm HI data available to understand the Galactic 
dust component of the maps. 



2. Observations 

We have undertaken a series of targeted observations of the Galactic HI component us- 
ing the 100-m Robert C. Byrd Green Bank Telescope (GBT) at NRAO. The GBT has 
an unblocked aperture with a high main beam efficiency and reduced far sidelobes com- 
pared with conventional telescopes. Accurate calibration of the GBT data is discussed 
in §0 
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Maps of Galactic HI in the 21cm line were made on-the-fly at constant Galactic 
latitude, b, or constant declination, <5, on a 3 '.5 grid, slightly finer than Nyquist sampling 
for the 9'1 beam of the GBT. The 2.5 MHz frequency-switching mode removes most 
of the continuum and filter response; a third order polynomial was fit to non-HVC 
channels outside of -100 km s" 1 and +25 km s" 1 to remove the remaining residual 
instrumental baseline. 

Some fields were observed several times to improve the sensitivity and allow us to 
estimate systematic effects brought about by stray radiation, spectral baselines, or other 
instrumental problems. Field centres, total integration times and resultant sensitivities 
are shown in Table Q] 



Table 1 . GBT high-latitude observations 



Field name 


Field centre 


Size 


tint 


Channel 


Scan Direction 




/ 


b 




(sec) 


noise (K) 




Polaris 


125? 00 


27? 50 


6° X 10° 


4 


0.18 


b 


KnotC 


138?99 


29=42 


2° x2° 


8 


0.12 


b 


Spider South 


135?70 


29? 80 


12° X 8.5° 


8/4 


0.11/0.17 


b 


NEP 


96? 40 


30?()0 


12° X 12° 


12 


0.09 


b 


KnotB 


147? 64 


34?77 


2° x2° 


8 


0.14 


b 


KnotN 


100?28 


35?50 


2° x2° 


4 


0.18 


b 


Ursa Major East 


155?80 


37? 00 


10.5° x 6° 


4 


0.16 


b 


Polaris North 


125? 00 


37? 40 


6° x 10° 


8/4 


0.12/0.16 


b 


draco 


92?24 


38?43 


5° x5° 


12 


0.11 


d 


Ursa Major 


143? 60 


38?5() 


9° x9° 


4 


0.15 


b 


KnotA 


152?62 


39?()0 


2° x2° 


8 


0.13 


b 


DRAOdeep 


134? 98 


40?()0 


10° x 10° 


8/4 


0.11/0.14 


b 


Nl 


85? 33 


44?28 


5° x5° 


8 


0.13 


b 


SP 


132?38 


47?50 


5° x5° 


8 


0.13 


b 


g86 


87? 94 


59?05 


5° x5° 


12 


0.11 


3 


bootesb 


55?98 


63? 87 


2.5° x4° 


12 


0.10 


S 


AG 


164? 85 


65?50 


5° x5° 


8 


0.13 


b 


bootesd 


56? 60 


65?61 


2.5° x4° 


12 


0.10 


6 


Necklace 


67? 82 


67?76 


2° x2° 


8 


0.12 


b 


bootes 


58? 11 


68?57 


4° x4° 


15/6 


0.10/0.16 


S 


bootesc 


60? 38 


71?45 


2.5° x4° 


12 


0.10 


S 


bootesa 


62?22 


73? 11 


2.5° x4° 


8/4 


0.16/0.17 




MC 


56? 84 


-81?50 


6° x5° 


8 


0.14 


b 



3. Data Reduction 

The GBT is not completely free from "stray" 21cm radiation. Although it lacks the 
scattering rings that arise from aperture blockage, it s till has a foreward spillover lobe 
( Lockman & Condonll2005l : iRobishaw & Heiles 12009) that contains several percent of 



the telescope's response at 21cm wavelength. In order to properly correct for the stray 
radiation, an empirical model of the telescope's sidelobes was developed. Properties of 
the main beam were established from electromagnetic calculations, and the sidelobes 
were measured from observations of the Sun (Figured]). The predicted stray radiation 
for a given observation, T stmy , that enters the receiver from directions away from the 
main beam, is determined by summing the amount of HI in every direction above the 
horizon (using the LAB HI survey) weighted by th e GBT response i n that direction. 
Complete information on the procedure is given in iBoothrovd etaD (l2010h : here we 
will summarize some of the results. 
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Figure 1 . GBT response pattern from spillover past the subrefiector, as measured 
from scans of the Sun. 



The antenna temperature scale was established through laboratory calibrations and 
observations of radio flux calibration sources. The two methods give results that differ 
by only a few percent. The GBT antenna temperature, T a , is corrected for atmospheric 
attenuation and converted to brightness temperature, Tb, adopting a main beam effi- 
ciency, rj m b - 0.88, determined from a theoretical calculation of the main beam profile. 
The final HI brightness temperature is then 



tray 



Vmb 



(1) 



where T* is the GBT antenna temperature corrected for atmospheric attenuation, and 
T stray is the calculated sidelobe contribution to the observed antenna temperature pro- 
file. 

Figure |2] shows the integrated HI from a field as observed over several different 
sessions with especially large and varying amounts of stray radiation (left), and the 
same data after the calibration, correction for stray radiation and baseline removal. The 
effects of the correction are evident in the reduced number of jumps in the intensity 
scale. 

3.1. General stability 

Robishaw & Heilesl ( 20091) developed an approximate model for the GBT beam and 



used this to correct their data for stray radiation. In so doing, they concluded that the 
GBT 21cm system was not stable, but had gain fluctuations on the order of 10-15%, 
which appeared to vary as a function of LST The more detailed model of the stray 
radiation we have developed allows us to test the gain stability with more confidence, 
and separate any true gain variability fro m other effects. The telescope, receiver, and 
IF systems were identical to those used by iRobishaw & Heilesl ((2009), but we used the 
GBT Spectrometer rather than the Spectral Processor as the detector. 

Over the course of the high Galactic latitude observations, which occured in many 
sessions over a three-year period, two fields were revisited regularly to test the con- 
sistency of the calibration: S6 and S8. After reducing these data using our correction 
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Galactic Longitude Galactic Longitude 

Figure 2. "Worst-case" low column density integrated GBT cube before stray and 
baseline removal (left) and after (right), showing a cleaner end product. 



for stray radiation, we find the measurements to be consistent to within a few percent, 
with no evidence of a correlation between HI brightness temperature and date or LST. 
The scatter in the measurements is entirely consistent with our estimates for the uncer- 
tainties arising from noise, spectral baselines, and residual stray r adiation. We do not 
reprod uce, nor can we understand, the "gain variations" reported by Robishaw & Heiles 
(2009). We also note that many tens of thousands of HI spectra have been mea- 
sured with the GBT by a number of observers, using both the Spe ctral Processor and 
Spectr o meter, without detecting significant gain variations (e.g., lLockman & Condon 
(120051) : IChvnoweth et all (l2008h ). 



3.2. A quick calibration check 



To confirm the accuracy of the final conversion to brightness temp erature (as in Equa- 
tion [J), we compare our mean integral of S6 and S8 with those of Williams! (1 19731) and 
IKalberla et all (1 19821) : £i 7^Av = 299 K km s" 1 , and 846 K km s _1 , respectively. The 
mean integrals of our GBT data are 289 K km s" 1 for S6, and 845 K km s -1 for S8 in 
remarkable agreement with both lWilliamsl (1 19731) and iKalberla et al.1 (1982). 



4. Results 



4.1. Data quality 

The integrated HI emission over the 12°by 12° NEP mosaic is shown in Figure [3] 
It is clear from this mosaic of observations at different LSTs (and thus very different 
positioning of the beam's sidelobes on the sky) that residual effects, which would appear 
as stripes or discordant spots, are small. 
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Figure 3. Total integrated HI intensity (zeroth moment) map of NEP cube 



4.2. Dust-gas correlation 

In the optically thin case the integrated line intensity, Whi (in K km s _1 ), can be con- 
verted to a column density, Nhi- 



Nhi = 0.018231V/// 



(2) 



where Nhi is in units of 10 20 cm" 2 . 

The GBT HI/IRIS 100 jum correlation is shown in Figure |4]for every independent 
observation, covering over 800 square degrees of high Galactic latitude sky. The cor- 
relation is quite strong and the slope (0.67 MJy sr" 1 ) is consistent with that found by 
Boulanger et al.l ([1996) (0.69 ± 0.03). However, there are clear deviations from the av- 



erage. At high column densities, an infrared excess can be attributed to the presence 
of mo lecular hydrogen (Boulan ger & Peraultlll988UJoncas et al.ll 19921 : iBoulange r et al. 
1996), but this does not explain the scatter seen throughout Figure 0] Does the scatter 
arise from real emissivity variations across the sky, or pos sibly emissiv i ty var iations 
along the line-of-sight? The former analysis is covered in iMartin et al.l (1201 Oh while 
the latter is discussed in the light of velocity tomography in the following section. 



4.3. Velocity tomography 

A number of the GBT fields were selected for the presence of IVC or HVC gas. If these 
clouds can be separated from the gas consistent with Galactic rotation (the 'local' gas) 
and from each other, a series of dust-emissivity-related coefficients {ax, b\, ...) can be 
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Figure 4. GBT Nhi correlation with IRIS 100/mi over more than 800 square de- 
grees of high Galactic latitude sky 

determined for each IR or sub-mm map, I A . An example of such a velocity separation 
is shown in Figure [5] The coefficients are then derived from solutions to 

h = a A + b A N h HI + c A N c m + (3) 
where component a is the cosmic infrared background. 




Galactic Longitude Galactic Longitude Galactic Longitude 



Figure 5. Local (left), IVC (centre) and HVC (right) components of the 21cm 
emission in the NEP field. Note differences in intensity scale. 

With knowledge of these emissivities for each given component (a,b,c,...) at a va- 
riety of wavelengths, a spectral energy distribution (SED) can be fit and dust properties 
for the given component can be determined. With the anticipated contribution from 
Planck in the sub-mm, this is on the verge of becoming an even more useful analysis 
of the dust content of, and dust evolution in ou r Galaxy. For now, this analysis has 



been done using IRIS (reprocessed IRAS) data (Miville-Deschene s & Lagachell2005l) 
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at 60 jim and 100 /urn. The sample field (NEP, as in Figure [5]) shows evidence of en- 
hanced 60 jum/lOO jum colour in the IVC as compared to the local gas as has been seen 
in the Draco IVC (Herbstmeier et al. fl993h - This is an indicator of dust processing with 
a slight increase in the relative contribution from small grains to the SED. In another 
field, Bootes, the IVC dust not only appears to be processed, but it appears to be ex- 
posed to a different (stronger) external radiation field. Understanding this in the context 
of the assumed constant interstellar radiation field at high Galactic latitude is key. Ap- 
parently, not all IVC gas is necessarily of the same origin, or at the least it cannot all be 
described as having an identical mixture of gas and dust, heated identically. 



5. The future 

In addition to studying the correlation between gas and dust, these GBT HI maps offer 
a large area over which to study the gas turbulence, the interaction between phases of 
the ISM and the properties of compact HVCs, amongst other things. This area of study 
promises to be quite fruitful over the next few years. 
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